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I*-  SPONSORING  MILITARY  ACTIVITY 

Aerospace  Medical  Research  Laboratory, 
Aerospace  Medical  Div. ,  Air  Force  Systems 
Cowand,  Wright-Patterson  AFB,  OH  45433 

Measurements  were  made  of  the  pilot  station  vibration  environment 
experienced  in  the  B-52  aircraft  during  low-altitude,  high-speed  flight. 

Data  were  obtained  of  the  longitudinal,  lateral  and  vertical  linear 
accelerations  and  the  roll  axis,  angular  velocities  and  accelerations.  The 
terrain  conto-.tr  following  flights  at  500-900  feet  above  the  flat  to  send- 
mountainous  regions  of  Louisiana  and  Arkansas  were  made  at  350  knots  true 
airspeed,  Data  analyses  included  probability  density  and  distribution  and 
the  auto-power  spectral  density  functions  in  addition  to  tests  for  station- 
arity,  randosness  and  normality.  The  individual  degree-of-freedom  results 
are  presented  in  graphical  and  tabular  form  and  in  general  indicate  that 
the  pitot  station  vibration  environment  produced  by  the  coafcination  of  gust 
response  and  maneuvering:  (1)  is  stationary  for  up  to  590  seconds;  (2)  is 
a  random  phenomena  that  does  not  satisfy  the  chi-square  goodness-of-fit  test 
for  Gaussian  distribution;  (3)  cannot  be  adequately  simulated  in  the  labora¬ 
tory  for  human  biomedical  tolerance  and/or  psycho-physiological  performance 
studies  using  only  vertical  axis  motion;  (4)  is  remarkably  in  agreement 
with  these  PSD's  calculated  for  the  lateral  and  vertical  axes  using  an 
aircraft  transfer  function  experimentally  derived  on  another  Drogram  and 
the  Dryden  gust  input  spectrum. 
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terrain  contour  flying  at  low-altitude,  high-speed  as  a  radar  defense 
penetiation  aid  increases  the  vibration  environment  experienced  in  an 
aircraft.  This  additional  stress  can  seri  ‘»sly  compromise  the  ability  of 
the  crew  members  to  satisfactorily  perform  their  assigned  duties.  Before 
human  biomedical  tolerance  and/or  psycho-physiological  performance  criteria 
can  be  established  by  using  laboratory  simulation  devices,  it  is  necessary 
to  first  define  tiic  vibration  environment  typical  of  these  military  type 
flight  operations.  Sincr  aircraft  design  characteristics  primarily  determine 
tiie  frequency  variations  in  the  vibration  environment  experienced  by  the  crew, 
the  aircraft  commonly  used  in  Air  Force  operations  have  been  categorized  into 
the  three  major  divisions  of:  (1)  small,  rigid  fighters;  (2)  large,  flexible 
bombers,  ( .’)  large,  rigid  bombers.  This  endeavor  concludes  the  study  under¬ 
taken  to  define  the  crew  compartment  vibration  environment  during  low-altitude, 
high-speed  flight  in  a  present  day,  large,  flexible  bomber. 

Appropriate  instruments  wore  installed  in  the  pilot  station  of  B-52F. 
model  aircraft  to  measure  and  record  the  vibration  environment  during  six 
low-alt i tude .  high -sneed  flights  over  flat  to  semi -mountainous  terrain. 
Measurements  were  made  of  the  linear  accelerations  in  the  longitudinal,  lateral, 
and  vertical  axes  and  the  roll  axis  angular  velocities  and  accelerations.  The 
acceleration  time  histories  were  analyzed  to  ascertain  their  amplitude  proba¬ 
bility  -nd  spectral  (frequency)  properties. 
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b'bcn  analyzing  random  vibration  data,  it  is  common  practice  to  assume 
the  data  i<;  statistically  stationary  and  exhibits  a  Gaussian  or  normal 
lmniitude  probability  distribution.  These  data  were  stationary  in  all  degrees- 
of- freedom  for  periods  of  12d-5yo  seconds,  but  did  not  meet  the  chi-square 
■toodness-of-fit  test  for  normality.  Over  the  freouency  range  (0.4  -  20  Hz) 
that  tip*  human  bodv  reacts  most  unfavorably,  power  spectral  density  analyses 
a  led  that  all  axes  contribute  significant  lv  to  the  overall  linear  and 
.'fi';;!  I -;»■  i;cc  ler:u  ion  environment  "resent  in  tnc  h-52  aircraft  during  low- 
'Ifi'  h!<  ,  hi  eh-spee.l  flight. 


•  f.*ined  accelerations  induced  by  maneuvers  greatly  increase  the 
>»  i  ,*e»>sit/  of  the  vibration  environment  and  are  probably  why  tiic  time 
;■  not  truly  random  amplitude  probability  density  functions, 

ip  i  mv  s i iv 1 1 n t i or  of  the  1..-V1S  vibration  environment  for  biomedical 

u-  •/”!  •  sycho-phvs iol pvi cn l  performance  studies  cannot  be  limited 
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INTRODUCTION 


The  increasing  i  rnortancr  being  placed  on  low-altitude,  high-speed  (IJVHS) 
missions  in  present  and  projected  aircraft  dictates  that  a  realistic  assess¬ 
ment  be  made  of  the  crew  station  vibration  environment  experienced  during 
such  flights.  Such  an  evaluation  is  necessary  in  order  that  valid  crew 
member  usycho-physiological  performance  and  tolerance  criteria  can  be 
established  through  the  use  of  in-laboratory  motion  simulation  testing. 

Since  basic  aircraft  design  characteristics  primarily  determine  the  vibration 
environment  felt  by  the  crow,  military  aircraft  have  been  categorized  into 
the  three  general  divisions  of:  (1)  small,  rigid  fighters;  (2)  large, 
flexible  bombers:  (3)  large,  rigid  bombers.  This  renort  presents  the  results 
of  measurements  and  analyses  made  of  the  crew  compartment  vibration  environ¬ 
ment  produced  by  the  response  of  a  large,  flexible  bomber  (B-52)  to  the 
combined  stresses  imposed  by  maneuvers  and  gusts  during  LAilS  flignt. 

An  in-flig'nt  measurement  program  was  performed  on  the  B-52  aircraft  in 
August  19o6 ,  during  which  l.(»  hours  of  LAJIS  data  was  obtained  on  six  flights 
conducted  at  350  knots  true  airsnecd  at  an  altitude  of  !>i>0-!h)'.>  feet  above 
slightly  mountainous  terrain,  i.incar  accelerations  in  the  longitudinal  f.\). 
lateral  (V)  and  vertical  (i.)  axes  and  angular  velocities  and  accelerations 
in  the  roll,  '’itch  and  yaw  axes  were  analog  recorded  as  amnlitude-time 
histories  for  subsequent  statistical  and  spectral  analvsis  in  each  axis, 
because  of  instrumentation  nroi: i eras  during  data  acquisition,  mo  Pitch  or 
yaw  axis  information  was  obtained  on  any  of  the  flights  and  nriy  flights 
1,  2,  4  and  (>  yielded  useable  data  in  the  other  degrees -of- freedom. 


SECTION  II 


MEASUREMENT  PROGRAM 

The  aircraft  used  in  this  endeavor  consisted  of  B-S2E  models  operationally 
assigned  to  the  Strategic  Air  Command  that  were  pressed  into  service  as 
test  bed  aircraft  for  a  large-scale,  low-level  flight  test  program  to  evaluate 
target  acquisition  capability.  The  low-altitude,  high-speed  vibration  data 
acquired  on  these  flights  were  obtained  as  a  "piggyback"  effort  to  the  primary 
mission  of  the  flight  test  program.  The  aircraft  were  typical  of  the  B-52 
design  with  the  exception  that  a  1  foot  diameter  and  6.7  foot  long  instru¬ 
mentation  pod  was  externally  centerline  mounted  to  the  underside  of  the 
aircraft  to  accomodate  the  test  equipment  required  for  the  primary  mission 
test  program.  Comparison  of  the  vibration  measurements  obtained  during 
these  tests  with  those  previously  made  on  another  B-52  program  indicates  that 
the  external  instrumentation  pod  did  not  perceptably  alter  the  crew  station 
vibration  environment  spectrum  in  either  the  lateral  or  vertical  axis.  The 
effect  of  the  pod  on  the  vibration  spectrum  in  the  other  degrees  of  freedom 
is  not  known  since  comparable  data  in  these  axes  were  not  available. 

The  vibration  transducer  assembly  used  in  this  program  was  rigidly 
mounted  to  the  cockpit  floor  in  front  of  the  instructor's  seat,  i.e., 
the  package  was  situated  between  and  slightly  behind  the  pilot  and  co-pilot 
seats.  Figure  1  is  a  cut-away  sketch  of  the  B-52  showing  the  locations 
within  the  aircraft  of  the  transducer  assembly,  tape  recorder  and  monitor/ 
control  console. 

INSTRUMENTAL  GN 

The  crew  station  six  degree  of  freedom  vibration  data  acquisition  system 
included  a  transducer  package,  individual  data  channel  visual  monitoring 


and  signal  conditioning  package  and  a  14  channel  analog  tape  recorder  with 
remote  control  capability.  The  transducer  assembly  contained  a  tri axial 
rate  gyro  device  and  three  linear  accelerometers  orthogonally  mounted  in 
the  configuration  shown  in  figure  2  so  as  to  align  the  respective  axes  with 
those  of  the  aircraft.  The  pitch  axis  rate  gyro  wa3  inoperative  throughout 
the  tests  and  the  yaw  axis  data  was  always  in  the  noise  floor,  therefore 
no  useful  data  was  obtained  in  these  degrees  of  freedom.  The  roll  rate 
gyro  measured  the  angular  velocity  in  the  independent  axis  over  the  ranges 
of  *100  degrees/sec.  To  supplement  static  calibrations  performed  with  an 
ordinary  rate  table,  the  gyros  were  dynamically  calibrated  at  the  Army 
Inertial  Guidance  and  Control  Laboratory  using  a  precision  oscillating  rate 
table.  The  linear  accelerations  were  measured  with  servo  type  accelerometers 
having  ranges  of  *lg  in  the  X  and  Y  axes  and  *Sg  in  the  Z  axis.  Calibration 
of  these  transducers  was  performed  using  a  laboratory  quality  low  frequency 
shake  table. 

To  optimize  the  quasi-static,  overall  data  signal  to  noise  ratio  that 
was  being  recorded  at  any  given  time,  the  analog  tape  recorder  input 
sensitivity  was  adjusted  to  a  very  low  level  voltage.  Each  data  channel 
included  a  40  dB  passive  attenuator  network  that  was  adjustable  in  10  dB 
increments.  A  trained  in-flight  operator  used  the  visual  monitoring  system 
to  determine  when  to  make  attenuator  setting  changes  end  thereby  assured  that 
the  maximum  overall  data  dynamic  range  was  attained.  In  this  manner, 
cocpensation  was  made  whenever  significant  changes  in  the  overall  vibration 
environment  was  roted  in  any  given  nxi».  Such  changes  in  the  overall  data 
signal  can  occur  during  L4HS  flight  over  different  types  of  terrain  during 
a  given  sortie.  However,  this  signal  conditioner  capability  did  not 


alleviate  the  signal  to  noise  ratio  limitations  imposed  by  the  tape  recorder 
on  the  instantaneous  dynamic  range  of  the  data  being  recorded. 

ROUTE  DESCRIPTION 

The  test  flights  for  this  study  were  flown  along  a  400  nautical  mile 
route  extending  from  England  Air  Force  Base  in  Louisiana  to  the  Ouchita 
Mountains  in  the  state  of  Arkansas.  The  portion  of  the  route  leading  into 
the  Ouchita  Mountains  is  rather  flat  with  the  highest  terrain  being  53S 
feet  above  sea  level  and  is  made  up  of  forests,  meadows,  rural  roads  and 
cultural  features  representative  of  small  to  medium  population  communities 
in  the  temperate  zone  worldwide.  That  part  of  the  route  within  the  Ouchita 
Mountains  is  primarily  tree  covered  and  rolling  type  terrain  with  the 
highest  peak  2,830  feet  above  sea  level.  For  navigational  purposes  the 
route  was  divided  into  an  East  Course  and  Nest  Course  with  no  discemable 
differences  in  the  general  topography  covered  by  the  two  segments  of  the 
route.  All  flights  were  performed  using  the  terrain  following  technique 
at  less  than  1,000  feet,  i.e.,  the  aircraft  were  flown  500-900  feet  above 
the  terrain.  The  test  pilots  were  instructed  to  maintain  as  best  they  could 
their  altitude  above  terrain,  airspeed  and  position  along  the  particular 
course  being  flown  at  that  time. 

TEST  FLIGHT  IDENTIFICATION 

On  each  flight,  data  was  recorded  for  a  maximum  of  sixteen  minutes  or  a 
distance  of  about  93  nautical  miles  along  the  route.  This  reord  time 
was  fixed  by  the  tape  capacity  and  recording  speed  of  the  recorder. 

Reference  1  states  that  the  properties  of  statistical  steticr.arity  and 
homogeneity  of  rather  mild  atmospheric  turbulence  might  be  expected  for 
large  scale  regions  of  up  to  85-90  nautical  miles,  but  for  cases  of  more 


severe  turbulence  as  can  be  encountered  over  mountainous  terrain  the 
distances  over  which  the  turbulence  exhibits  stationarity  are  usually 
somewhat  smaller.  Since  most  of  the  data  acquired  on  these  flights  was 
over  semi -mountainous  terrain,  before  any  spectral  or  amplitude  probability 
analyses  can  be  performed  on  the  data  it  was  necessary  to  sub-divide  each 
flight  record  into  segments  that  satisfied  the  statistical  property  of 
stationarity. 

Reference  2  presents  in  rigorous  terms  the  general  concept  of  stationarity 
and  the  verification  of  time  invariant  statistical  properties  theoretically 
required  for  proof  of  stationarity.  The  reference  shows  such  proof  of 
stationarity  is  clearly  not  feasible  in  practical  term3  since  it  would  require 
an  infinite  number  of  probable  statistics.  A  practical  test  for  stationarity 
is  developed  in  the  text  based  on  the  following  assumptions:  (1)  proof  of 
self  stationarity  for  individual  sample  records  can  be  accepted  as  proof 
of  stationarity  for  the  random  process  from  which  the  time  histories  were 
obtained;  (2)  verification  of  weak  stationarity,  that  is  the  mean  value 
and  autocorrelation  function  are  not  time  dependent,  will  be  acceptable  for 
most  practical  types  of  analyses  and  applications;  (3)  the  sample  record 
lengths  are  very  long  compared  to  the  random  fluctuations  of  the  data; 

(4)  proof  of  stationarity  of  the  mean  square  value  can  be  accepted  as 
proof  of  stationarity  of  the  autocorrelation  function.  This  last  assumption 
is  usually  valid  since  the  mean  square  value  is  equal  to  the  autocorrelation 
function  at  X  •  o  and  it  is  highly  unlikely  that  non -stationary  data  will 
have  a  time-varying  autocorrelar ion  function  that  decs  not  also  vary  at 
*C  •  o.  The  test  for  stationarity  consisted  of  merely  applying  the  reverse 
arrangement  or  trend  test  to  the  mean  and  mean  squaie  values,  obtained  for 


negjuo! 


each  consecutive  ten  second  interval  of  the  acceleration  tine  history. 

Each  flight  was  first  sub-divided  whenever  the  in-flight  operator  made 
an  attenuator  setting  change  in  the  signal  conditioning  system  during 
data  acquisition.  Each  flight  was  then  further  subdivided  when  the  test 
for  stationarity  failed  in  any  axis  at  the  0.05  level  of  significance. 
Table  I  lists  the  flight  segment  identification,  aircraft  number,  course 
designation  and  total  time  in  seconds  that  the  data  satisfied  the  test  for 
stationarity  in  ail  degrees  of  freedom.  Figures  3-11  show  the  terrain  and 
flight  profiles  above  mean  sea  level. 


TABLE  I 


FLIGHT  IDENTIFICATION 


IDENTIFICATION 

NUMBER 

AIRCRAFT’ 

NUMBER 

COURSE 

DESIGNATION 

STATIONARY 
PERIOD  (SEC) 

1-11 

640 

WEST 

230 

1-1-2 

640 

WEST 

130 

1-1-3 

640 

WEST 

130 

1-2-1 

640 

WEST 

140 

1-2-2 

640 

WEST 

140 

1-2-3 

640 

WEST 

120 

2-1 

640 

EAST 

120 

2-2 

640 

EAST 

230 

2-3-1 

640 

EAST 

140 

2  * 

640 

EAST 

330 

4-2 

127 

WEST 

280 

4-3 

127 

WEST 

380 

4-4 

127 

WEST 

170 

6-2 

EAST 

S90 

6-3 


634 


EAST 


230 


SECTION  III 


DArA  ANALYSES 

To  completely  define  any  multiple  degree  of  freedom  vibration  environment, 
it  is  necessary  to  determine  the  amplitude,  frequency  and  phase  characteristics 
of  the  data.  Since  the  immediate  application  of  this  data  is  to  use  it  as 
input  for  in-laboratory  simulation  of  the  environmen  psycho-physiological 

human  tolerance  testing  and  the  SIXMOUF,  motion  simulation  device  that  will 
be  used  in  such  studies  currently  has  an  independent  aegree  of  freedom 
control  capability,  amplitude-time  and  frequency  analyses  were  performed  in 
each  axis.  Joint  probability  and  cross-spectral  analyses  to  obtain  phase 
information  were  not  warranted  at  this  time. 

PROBABILITY  DENSITY  AND  DISTRIBUTION 

Once  it  has  been  determined  that  the  data  is  at  least  weakly  self 
stationary,  then  a  test  to  ascertain  randomness  is  in  order.  Such  tests 
usually  involve  analysis  techniques  that  assume  the  data  are  random.  The 
presence  of  periodic  or  near-periodic  components  can  be  detected  by  visual 
inspection  of  the  autocorrelation,  power  spectral  density  and/or  the  amplitude 
probability  density  functions.  Although  not  as  powerful  as  the  other  two 
methods,  use  of  the  amplitude  probability  density  function  was  chosen  fer 
the  following  reasons:  (1)  the  autocorrelation  function  docs  not  yield  any 
additional  information  over  that  obtained  from  the  power  spectral  density 
function  since  one  is  merely  the  Fourier  transform  or  inverse  transform 
of  the  other;  (2)  power  spectral  density  analysis  will  be  performed  on  the 
data  any  way  since  a  measure  of  the  frequency  composition  must  be  made  to 
supplement  any  amplitude-time  analysis  performed;  (3)  an  analysis  of  the 
amplitude-time  distribution  is  of  great  value  in  itself  in  the  description 


8 


of  a  given  vibration  environment. 


Digital  techniques  were  used  to  compute  the  normalized  amplitude 
probability  density  and  cumulative  distribution  functions  in  the  X,  Y.  Z 
and  roll  axes  for  each  flight  segment  that  exhibited  stationarity.  Each 
acceleration  time  history  was  digitized  at  a  sampling  rate  of  loo  points 
per  second  and  the  number  of  points,  N.,  that  fell  within  each  of  60  equal 
amplitude.  Ay,  intervals  was  computed.  Knowing  the  total  number  of  data 
points,  N,  in  the  sample  time  history,  the  amplitude  prooability  density 
function  is  defined  as: 

N; 

Probability  Density  Function  »  (y)  «  w^y  (2 

Ihe  probability  density  function  defines  the  probability  that  the  instanta¬ 
neous  value  lies  within  the  range  of  y  and  (y  ♦  4  y) .  the  probability  that 
the  instantaneous  value  is  less  than  or  equal  to  some  specified  value  y  is 
the  cumulative  probability  distribution  function,  and  it  is  defined  as: 

A 

Probability  Distribution  Function  »  P^  (y)  «  (y)  (2J 

The  computer  program  used  in  this  study  computes  the  normalized  probability 
density  and  distribution  functions  by  dividing  the  instantaneous  values  by 
the  standard  deviation.  *  The  standard  deviation  is  the  rms  variotion 
about  the  a»an  value  and  can  be  found  from  the  variance  which  is  defined  by: 

Variance  ■  (T*  *  (rms  value)2  -  (mean  value)2  (3\ 

Since  the  computer  program  forces  the  mean  value  to  zero,  then  it  is  seen 
from  equation  3  that  the  standard  deviation  is  equal  to  the  rms  value  of 
the  tine  history.  Figure  12  shows  the  scaling  used  in  all  of  the  probability 
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density  and  distribution  computations  performed  in  this  study.  Figures  13-72 

are  the  normalized  probability  density  and  distribution  functions  of  the 

X,  Y,  Z  and  roll  axis  accelerations  for  each  flight.  The  legend  shown  below 

each  graph  indicates  the  number  of  equal  amplitude  intervals  used  in  the 

computations  and  the  voltages  corresponding  to  the  root-mean-square  and 

mean  values.  Table  II  gives  for  each  flight-axis  the  root-mean-square 

2 

acceleration  value  in  g  and  radian/(sec)  units  after  taking  into  account  all 
transducer  sensitivities  and  data  acquisition  and  reduction  correction 
factors. 

POWER  SPECTRAL  DENSITY 

Because  of  the  biomedical  considerations  involved  in  the  expected 
application  of  the  I.AHS  data  collected  during  this  program,  the  frequency 
composition  of  the  data  is  of  utmost  importance.  Recent  improvements  have 
been  made  in  the  stability  and  reliability  of  analog  devices  designed  to 
measure  the  power  and  cross  power  spectral  density  functions.  Commercially 
available  analog  systems  can  now  be  used  to  perform  such  analyses  as 
accurately  as  v.as  previously  available  only  by  digital  means.  Historically 
speaking,  swept  spectrum  analyses  have  received  widespread  acceptance  for 
the  frequency  analysis  of  random,  neriodic,  or  mixed  data  samples  that  are 
either  stationary  or  transient  in  nature.  Their  appeal  has  been  mainly 
due  to  their  versatility  in  being  able  to  properly  analyze  a  wide  variety 
of  data  signals,  especially  acoustic  and  vibration  environments  at  the 
higher  frequencies.  Numerous  studies  of  the  effects  of  sinusoidal  and 
random  vibrations  on  man  have  shown  the  human  body  reacts  most  unfavorably 
to  those  accelerations  at  frequencies  less  than  20  Hz.  Primarily  due  to 
advancements  in  filter  and  detector  design,  frequency  analyses  can  now  be 
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FLIGHT 

IDENTIFICATION 

1-1-1 

1-1-2 

1-1-3 

1-2-1 

1-2-2 

1- 2-3 

2- 1 
2-2 
2-3-1 
2-3-2 
4-2 
4-3 
4-4 
6-2 
6-3 


CREW 

TABLE  II 

STATION  RMS  ACCELERATION  VALUES 

g  UNITS 

RAD/ SEC 2  UNITS 

X  AXIS 

Y  AXIS 

Z  AXIS  | 

|  ROLL  AXIS  | 

.03 

.06 

.IS 

.5 

.02 

.06 

.16 

.7 

.02 

.06 

.17 

.6 

.03 

.OS 

.14 

.5 

.04 

.OS 

.13 

.6 

.02 

.04 

.12 

.6 

.02 

.06 

.17 

1.5 

.02 

.06 

.18 

1.5 

.02 

.05 

.16 

1.4 

.02 

.05 

.16 

1.5 

.02 

.06 

.16 

.6 

.OS 

.13 

.18 

.6 

.05 

.07 

.22 

.7 

.03 

.05 

.15 

.6 

.07 

.05 

.16 

.6 

performed  with  swept  spectrum  analyzers  at  the  low  frequencies  where  the 
limiting  factors  in  the  analyses  are  certain,  inviolate,  underlying 
statistical  principles  and  not  the  performance  characteristics  of  the 
analyzer. 

There  is  an  error  associated  with  any  measurement  or  analysis  effort 
dealing  with  random  data  since  the  data  samples  have  a  finite  time  duration 
and  can  only  be  described  by  statistical  properties.  In  addition,  most  of 
these  descriptive  quantities  arc  merely  estimates  or  probabilities  of 
occurrence  rather  than  exact  values.  In  general,  the  statistical  number 
of  dcgrees-of-freedom  is  a  measure  of  the  accuracy  of  the  results  obtained 
using  a  given  frequency  analysis  procedure  and  is  defined  by: 

n  =  2  BT  (4) 

where  B  is  the  effective  bandwidth  of  the  filter  employed  and  T  is  the 
sample  record  length  or  integration  time.  The  lengths  of  time  that  all  axes 
exhibited  stationarity  as  listed  in  Table  I  were  used  as  the  value  of  T  for 
each  flight  of  the  R-52.  From  equation  4,  once  T  has  been  established,  some 
trade-off  must  be  made  between  the  desired  filter  bandwidths  and  the  accuracy 
obtained  for  the  analyses.  To  ascertain  all  major  resonance  and  null  regions 
present  in  the  B-52  acceleration  data  were  being  defined,  one  data  sample 
was  analyzed  using  a  0.05  llz  filter.  As  a  reasonable  compromise  between 
frequency  resolution  and  measurement  error,  all  of  the  power  spectral 
density  measurements  were  made  using  a  0.2  Hz  filter.  This  means  even  with 
perfect  instrumentation  and  careful  calibration,  there  will  be  an  uncertainty 
associated  wi ih  these  mower  spectral  density  results  purely  because  ot  the 
statistical  limitations  involved  when  dealing  with  random  data.  Using 
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standard  statistical  procedures  it  is  possible  to  specify  the  confidence  level 
or  percent  probability  a  given  measured  value  does  lie  within  some  specified 
accuracy  limits.  For  example,  considering  flight  1-1-1  at  the  00%  confidence 
level  with  n  =  2  BT  =  2  X  0.2  X  230  =  92  statistical  degrees  of  freedom, 
one  would  find  00%  of  the  time  the  measured  power  spectral  density  is  within 
♦23.0%  or  _+1.8  dB  of  the  actual  value.  In  addition  to  this  inaccuracy  due 
to  the  limited  number  of  degrees  of  freedom,  one  must  also  account  for 
various  data  acquisition  and  reduction  uncertainties  nresent  in  the  PSD 
analyses.  Due  to  various  instabilities  in  the  total  data  acquisition  system, 
these  recorded  acceleration  time  histories  inherently  have  a  *S%  error 
associated  with  them.  An  additional  0.5%  uncertainty  in  the  PSD  results 
stems  from  the  accuracy  limitations  of  the  power  detector  used  in  the 
swept  spectrum  analyzer.  Table  III  lists  the  total  estimated  accuracy 
limits  at  the  90%  confidence  level  that  can  be  expected  of  the  PSD  curves 
developed  in  this  study. 

The  PSD  curves  shown  in  figures  73-87  were  obtained  with  a  fixed  filter 
bandwidth  of  0.2  Hz  and  time  integration  being  applied  to  the  power  detector 
output  over  the  entire  sample  record  length. 

Figures  88-92  show  for  each  individual  axis  all  of  the  crew  station 
power  spectral  densities  that  were  measured.  The  spread  in  the  data  shown 
in  figures  88-92  was  expected  since  it  includes  not  only  the  variations 
due  to  gust  input,  nut  also  those  due  to  terrain,  aircraft  number,  maneuvering, 
etc.  For  example,  figures  93-97  show  the  variability  in  the  power  speUrai 
density  in  each  axis  for  various  segments  of  Flight  1. 


TABLE  III 


PSD  ANALYSIS-MEASUREMENT  ERROR  AT  90*  CONF! HENCE  LEVEL 


FLIGHT  NO. 

DEGREES  OF  FREEDOM 

ERROR 

(IB  ERROR 

1-1-1 

92 

31.1 

2.4 

1-1-2 

52 

39.1 

2.9 

1-1-3 

52 

39.1 

2.9 

1-2-1 

36 

37.9 

2.8 

1-2-2 

56 

37.9 

2.8 

1-2-3 

48 

40.4 

3.0 

2-1 

48 

40.4 

3.0 

2-2 

92 

31.1 

2.4 

2-3-1 

56 

37.9 

2.8 

2-3-2 

132 

27.1 

2.2 

4-? 

112 

28.8 

2.2 

4-3 

152 

25.7 

2.0 

4-4 

98 

30.3 

2.3 

6-2 

236 

22.0 

1.7 

6-3 

92 

31.1 

2.4 

14 
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SECTION  IV 


COMPARISON  OF  EXPERIMENTAL -HYPOTHETICAL  POWER  SPECTRAL  DENSITIES 


The  fundamental  relationship  between  the  innut  and  response  power 
spectral  density  for  any  constant  parameter  linear  system  is  given  by: 

PSD  =  M2  X  PSD.  .  (S) 

where  JtJ2  is  the  gain  factor  or  transfer  function  of  the  system  or  "black 
box"  in  question  and  is  defined  as  ^TR*^  ♦  T^  wl  ere  TR  and  Tj  are  the 
real  and  imaginary  components  of  the  complex  frequency  response  function. 

From  reference  3,  experimentally  derived  Y  and  Z  axis  crew  station 
transfer  functions  were  extracted  for  a  B-52  aircraft  flying  under  test 
conditions  very  similar  to  those  present  during  this  data  acquisition  program. 
The  aircraft  gross  weight  and  airspeed  were  nearly  identical,  but  the  flight 
altitude  differed  considerably.  However,  the  flight  altitude  should  primarily 
be  reflected  in  the  level  of  turbulence  and  response  seen  at  the  crew  station, 
but  it  should  not  seriously  affect  the  spectral  characteristics  of  the  transfer 
function. 


To  obtain  a  hypothetical  power  spectral  density  response  function,  the 
values  for  the  crew  station  transfer  function  were  multiplied  on  a  point  by 
point  basis  by  the  power  spectral  density  of  the  gust  input  spectrum,  which 
was  empirically  derived  by  Dryden  as: 


.  <r*L  '  + 

gust  input  Tf  v  7  ~  x2-]A 

'  w 
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wstore  C T  is  tae  root  '■’can  square  velocity  of  the  gust  inrmt  in  ft/sec,  u)  is  the 
frequency  in  radians/sec ,  1.  is  the  turbulence  scale  factor  in  feet  and  V  is 
the  true  airspeed  of  the  aircraft  in  ft/sec  units.  Reference  3  states  that 
tno  same  gust  input  expression  can  lie  usee  tc  realistically  represent  hotli 
the  lateral  and  vertical  turbulence  spectra.  This  means  that  using  L  =  SUO 
feet  anu  V  =  r.91  l't/scc.  equation  6  can  be  multiplied  at  discrete  frequencies 
by  the  annronriate  transfer  function  value  to  obtain  both  the  Y  and  Z  axis 
normalized  crew  station  response  nower  spectral  density.  Figures  98  and  99 
snow  t.ic  Y  ami  Z  axis  nower  spectral  densities  measured  on  Flight  1-1-1  anu 
those  calculated  using  tiie  experimental  transfer  functions  from  reference  3 
and  the  hrvden  gust  innut  snectrun  with  £T“  =  19  (ft/sec)“  for  the  lateral 

>  7 

response  ami  y  -  33  (ft/sec}“  for  the  vertical  response. 
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SECTION  V 


CONCLUSIONS 

As  shown  in  Table  I,  the  vibration  environment  in  the  X,  Y,  Z  and  Roll 
axes  exhibited  the  statistical  property  of  stationarity  for  periods  of 
120-590  seconds.  The  distances  flown  during  these  neriods  of  data  stationarity 
varied  from  13. £  to  63  miles.  Reference  4  showed  that  an  F-4C  aircraft  was 
able  to  travel  20-75  miles  during  periods  of  vibration  data  stationarity 
while  also  terrain  contour  flying  over  semi -mountainous  regions  but  at  much 
faster  speeds  than  the  B-52  used  in  this  study.  These  distances  flown  during 
data  stationarity  are  in  excellent  agreement  with  each  other  and  are  consist¬ 
ent  with  the  maximum  of  100  miles  predicted  in  reference  1  for  large  scale 
regions  where  the  properties  of  stationarity  and  homogeneity  of  turbulence 
might  be  expected  to  aoply  during  non-storm  conditions.  They  are  quite  less 
than  the  approximately  130  mile  distances  covered  during  data  stationarity 
as  mentioned  in  reference  3.  However,  special  consideration  should  be  given 
to  the  fact  that  during  the  B-52  flights  conducted  in  reference  3,  only 
light  turbulence  was  encountered  and  the  flights  were  made  at  500-1000  feet 
above  the  highest  terrain  obstacle  rather  than  terrain  contour  flying  as 
was  done  in  this  study.  Detailed  examination  of  our  B-52  acceleration  time 
histories  and  route  profiles  revealed  that  the  periods  of  data  stationarity 
were  generally  terminated  when  the  aircraft  underwent  a  sustained  maneuver 
or  flew  over  terrain  where  the  general  characteristics  changed  considerably, 
e.g..  from  relatively  flat  terrain  to  semi -mountainous .  Therefore,  it  is 
not  surprising  that  the  periods  of  stationarity  in  this  study  may  be  less 
than  those  experienced  in  other  programs  after  it  is  recognized  that  aircraft 
maneuvers  are  prevalent  during  low-altitude,  high-speed  flight  and  that 
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regions  of  homogeneous  turbulence  are  general! y  limited  to  covering  terrain 
of  similar  characteristics. 

Visual  inspection  of  the  probability  density  functions  in  figures  13-72 

[  indicates  the  vibration  time  histories  in  any  one  axis  are  generally  random 

[ 

in  nature,  but  the  vast  majority  of  them  contain  one  or  more  high  probability 
components  at  specific  amplitude  intervals.  The  probability  density  function 
is  a  plot  of  the  percentage  of  the  total  time  spent  at  a  given  .amplitude 
level.  Obviously,  maneuvers  would  tend  to  produce  small  amplitude  changes 
over  a  long  period  of  time  compared  to  the  length  of  time  spent  at  a  given 
level  due  to  purely  random  fluctuations.  The  particular  amplitude  interval 
containing  these  accelerations  produced  by  a  maneuver  would  vary  according 
r  to  the  maneuver  characteristics,  and  would  explain  why  the  spikes  seen  in 

these  probability  density  functions  occur  at  various  amplitude  intervals. 

\ 

j  Since  maneuvers  while  terrain  contour  fiying  are  not  performed  periodically 

i  i 

and  the  overall  shapes  of  these  probability  density  figures  are  not  similar 

t  to  those  associated  with  mixed  random  and  periodic  phenomena,  it  is  hypothe¬ 

sized  that  the  extremely  high  and  low  amplitude  intervals  present  in  these 
data  are  the  result  of  aircraft  maneuvers. 

j  j 

i; 

When  dealing  with  turbulence  associated  studies,  it  is  common  practice 

i 

to  assume  the  measured  data  exhibits  a  normal  or  Gaussian  distribution. 

i 

This  is  because  the  Central  Limit  Theorem  suggests  that  any  phenomenon  that 
is  the  result  of  numerous  random  effects  will  tend  to  have  &  Gaussian 

j 


distribution  itself.  However,  strong  nonlinear  operations  can  produce 
deviations  from  the  ideal  Gaussian  form.  Negative  results  were  obtained 
applying  the  chi-square  goodness -of- fit  test  for  normality  to  the  five 


probability  density  functions  that  appeared  to  most  closely  follow  the 
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hell-shaped  density  function  characteristic  of  Gaussian  distributed  phenomena. 
Based  on  these  findings,  the  B-52  vibration  environment  experienced  at  the 
crew  station  during  Jow-altitude,  high-speed  flight  cannot  be  considered 
Gaussian  in  accordance  with  the  statistical  definition  of  the  term. 

Fro*  the  viewpoint  of  in-laboratory  simulation  of  the  crew  compartment 
low-altitude,  high-speed  in-flight  vibration  environment,  the  most  significant 
conclusion  drawn  from  this  study  can  be  readily  derived  from  figures  73-87 
showing  the  power  spectral  density  in  each  individual  axis.  Most  human 
biomedical  tolerance  and  psychophysiological  performance  studies  concerning 
the  effects  of  such  complex  aircraft  vibration  environments  that  have  been 
conducted  in  the  past  were  limited  to  consideration  of  the  vertical  axis 
accelerations  only.  However,  figures  73-87  unquest ion ably  show  that  above 
3  Hz  the  overall  linear  vibration  environment  measured  in  the  crew  compart¬ 
ment  was  not  solely  dependent  upon  the  vertical  axis  accelerations,  but  that 
the  lateral  axis  levels  were  equal  to  or  greater  than  those  in  the  2  axis. 

In  addition,  above  approximately  10  Hz,  the  longitudinal  axis  levels  were 
nearly  equivalent  to  those  seen  in  the  lateral  and  vertical  axes.  As  reported 
in  reference  4,  it  was  significant  to  ascertain  that  the  linear  accelerations 
in  all  axes  must  be  considered  during  any  research  programs  on  human  tolerance 
and/or  performance  criteria  and  the  in-laboratory  simulation  of  these  complex 
vibration  environments. 

F.xcludi-.g  the  very  low  frequency  or  phugoid  vibrations,  the  following 
characterizations  concerning  the  spectral  content  of  the  acceleration 
environment  in  each  axis  of  the  B-52  can  he  made  from  figures  88-92: 

X  axis  -  null  regions  located  at  1.2,  1.8,  2.5  and  3.6  Hz  separated  by 
rather  sharp  resonances  at  1.6,  2.0,  3.0  and  5.0  Hz  followed  by  a  mild  10  dB 


19 


decay  to  15  Hz  and  then  increasing  to  20  Hz ;  Y  axis  -  primary  peak  at  3.6  Hz 
with  other  resonances  located  at  1.3,  2.S  and  4.7  llz  followed  by  a  30  dB 
decay  out  to  20  llz.  These  resonances  are  separated  by  narrow  null  regions 
located  at  O.S.  2.0,  2.8  and  4.2  llz:  Z  axis  -  null  regions  located  at  1.2, 

2.3,  3.5,  4.3  and  6.0  Hz  separated  by  resonances  at  1.9,  3.0,  3.8,  5.7  and 
6.2  llz  followed  by  a  30  dB  decay  out  to  20  llz;  Roll  axis  -  null  regions 
located  at  1.0,  2.3,  3.3,  4.2  and  6.0  llz  separated  by  narrow  resonances  at 

1.3,  2.8,  3.6  and  5.0  Hz  with  a  broad  resonance  region  located  about  9.0  Hz 
followed  by  a  mild  8  dB  decay  to  12  llz  and  then  deteriorating  into  a  broadband 
signal  out  to  20  llz  with  no  discemabie  null  or  resonance  regions.  As 
expected,  the  roll  axis  angular  velocity  and  acceleration  environments 
exhibit  identical  resonant  and  nodal  characteristics.  Comparison  of  the 

roll  axis  spectra  does  show  the  angular  acceleration  levels  increase  with 
increasing  frequency.  This  increase  in  amplitude  with  increasing  frequency 
is  inherent  in  the  differentiation  process. 

As  seen  in  figures  88-90,  the  linear  acceleration  PSD  in  the  longitudinal, 
lateral  or  vertical  axis  only  varied  by  about  5  dB  for  all  flights.  In 
view  of  the  widely  different  terrain,  gust  input  and  pilot  maneuver  conditions 
under  which  the  data  was  acquired,  this  spread  in  the  PSD  spectra  was 
surprisingly  small. 

Figures  95-97  point  out  the  even  more  astounding  conclusion  that  the 
variability  in  the  PSD  spectra  within  a  given  flight  is  equal  to  the 
variations  seen  when  comparing  different  aircraft  flying  over  different 
terrains  on  different  days.  These  data  show  the  B-52  linear  acceleration 
crew  station  environment  experienced  on  one  flight  is  probably  no  worse 
than  S  dB  different  than  that  experienced  or.  any  other  flight  made  under  the 


same  general  flight  profile  conditions. 

Figures  98  and  ‘>9  show  the  almost  exact  spectral  agreement  between  our 
Flight  1-1-1  lateral  and  vertical  PSD  curves  and  those  calculated  using  the 
experimental  transfer  functions  from  reference  3  and  the  empirically  derived 
Dryden  gust  input  spectrum.  Since  no  measurement  of  the  gust  input  intensity 
was  made  on  any  of  the  flights  in  this  study,  the  values  chosen  for  (T  in 
figures  98  and  99  were  selected  to  yield  the  best  possible  amplitude 
agreement  with  the  Flight  1-1-1  curves.  However,  the  (T values  used  are 
quite  consistent  with  those  expected  during  low-altitude,  high-speed  flight 
over  semi-mountainous  terrain.  In  addition,  the  (f value  only  controls  the 
relative  position  of  the  PSL)  curve  along  the  amplitude  scale  and  in  no  way 
does  the  choice  of  the  <T value  affect  the  spectral  characteristics  of  the 
PSD  response  function. 

In  summary,  the  more  important  conclusions  drawn  from  this  effort  were: 

(1)  The  time  periods  over  which  all  axes  exhibited  the  statistical 
property  of  sel f-stationarity  were  limited  to  120-590  seconds  due  to 
sustained  aircraft  maneuvers  or  large  changes  in  the  intensity  of  the  gust 
input  spectrum  from  flying  over  terrain  with  quite  varied  general  landscape 
characteristics . 

(2)  Sustained  accelerations  induced  by  maneuvers  strongly  contribute 
to  the  overall  intensity  of  the  vibration  environment  and  are  probably  why 
the  time  histories  do  not  exhibit  truly  random  probability  density  functions. 

(3)  The  common  assumption  of  Gaussian  distribution  docs  not  seem 
warranted  for  these  data  since  all  cases  evaluated  by  the  chi-square 
goodness -of- fi t  test  for  normality  failed  at  the  0.93  level  of  significance. 

(4)  All  axes  contribute  significantly  to  the  total  vibration  environment 

:i 


experienced  by  the  crew  of  B-52  aircraft  flying  low-altitude,  high-speed 
missions. 

(5)  In- laboratory  simulation  of  the  LAHS  vibration  environment  for 
biomedical  tolerance  and/or  psycho-physiological  ncrformance  studies  cannot 
be  limited  to  motion  in  only  one  degree-of- freedom. 

(6)  variations  in  aircraft  or  terrain  did  not  affect  the  pilot  station 
PSD  in  the  linear  acceleration  degrees -of- freedom  by  more  than  5  dB  for 
frequencies  up  to  20  Hz. 

(7)  These  measured  lateral  and  vertical  axis  PSD  response  functions 
are  nearly  identical  in  spectral  content  to  those  calculated  using  aircraft 
transfer  functions  experimentally  derived  under  different  flight  conditions 
on  another  program  and  theoretical  gust  input  spectra. 
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Figure  3.  Flights  1-1-1,  1-1-2  and  1-1-3  Profiles 
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'igure  5.  Flights  2-1  and  2-2  Profiles 
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Figure  6.  Flights  2-3-1  and  2-3-2  Profiles 
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Flight  4-2  Profile 
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Figure  10.  Flight  6-2  Profile 
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Figure  12.  Probability  Scaling 
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!■  injure  iia.  Flight  1-i-I  X  Axis  Probability  Density 
(See  figure  12  for  sealing) 


Figure  16=1.  Fligh*  i  —1—1  SL.  II  Aucelerat.'i.-n  Irbability 


I6h.  Flight.  1-1-1  Roii  Acceleration  Probability  Distribution 
(See  Figure  12  tor  scaling) 
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Figure  18b.  Flight  1-1-2  ¥  Axis  Probability 
(See  figure  12  for  scaling) 
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Figure  20b.  Flight  1-1-2  Roll  Acceleration  Probability  Distribution 
{See  figure  12  .'.'or  scaling) 
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Figure  2ib.  Flight  1-1-3  X  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 


gure  22a.  Plight  i-l-j  Y  Axis  Probability  Densi 
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Figure  23b.  Flight  1—1 —3  Z  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 
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Figure  24b.  Plight 1-1-3  Roll  Acceleration  Probability  Distribution 


‘igurc  25a.  Flight  1-2-1  X  Axis  Probability  Density 
(oee  figure  12  for  scaling) 
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Figure  25b-  Flight  1-2-1  X  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 


'igure  26a.  Flight  1-2-1  Y  Axis  Probability  Density 
(t:ee  figure  12  for  scaling) 


Figure  26b.  Flight  1-2-1  Y  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 


2?a.  Flight  1-2-1  Z  Axis  Probability  Density 
(bee  figure  12  for  scaling) 


Figure  27b.  Flight  1-2-1  Z  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 
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Figure  28b.  Flight  1-2-1  Rcll  Acceleration  Probability  Distribution 
(See  figure  12  for  scaling) 


Figure  29a.  Flight  1-2-2  X  Axis  Probability  Density 
(i>e  figure  12  for  scaling) 


Figure  29b.  Flight  1-2-2  k  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 


Figure  3Ca.  Flight  1-2-2  Y  Axis  Prob. 

(See  figure  12  for  scalii 


Figure  30b.  Flight  1-2-2  Y  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 
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Figure  31b.  Flight  1-2-2  Z  A>J.i  Probability  Distribution 
(i?ee  figure  12  for  scaling) 
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Flight  1-2-2  Roll  Acceleration  Probability  Density 
(Dee  figure  12  for  scaling) 


p5^f^32b^FWgh^^2^2RoTI  Acceleration  Probability  Distribution 
(Sec  figure  12  for  scaling) 


igure 


pq  i  F  o‘oo  • 


INMIHRMPPJ* 


WMPIH  M  4  i IIPU* Ji P,  ft 


Figure  34b.  Flight  1-2-3  Y  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 


Figure  35a.  Flight  1-2-3  -  -V  is  Probability  Density 
(See  figure  12  for  scaling) 


Figure  35b.  Flight  1-2-3  1  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 
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gure  36b.  Flight  1-2-3  Roll  Acceleration  Probability  Distribution 
(See  figure  12  for  scaling) 


Figure  37h.  Flight,  2-1  X  Axis  IroUbility  Density 
(liee  figure  12  for  scaling) 
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Figure  37b.  Flight  i:  •'  X  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 


figure  3^a.  Plight  2-1  Y  Axis  i  robatdlit.y  Density 
(i'«e  figure  12  for  scaling) 
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'igure  38b.  Flight  2-1  Y  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 
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Figure  40a.  Flight  2-1  Roll  Accel ei-ation  Probability  Der. 
(w'ee  figure  12  for  scaling) 
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Figure  40.  Flight  2-1  Roll  Acceleration  Probability  Distribution 
(See  figure  12  for  scaling) 


Figure  4.1a.  Flight  2-2  X  Axis  Probability  Density 
(lee  figure  12  for  scaling) 
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Figure  41b.  Flight.  2-2  X  Axis  Probability  Distribution 
"igure  12  for  sealing) 
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re  u<b.  Flight  2-2  Z  Axis  hrobability  Distribution 
(See  figure  12  for  scaling) 
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Figure  U*b.  Flight  2-2  Roll  Acceleration  Probability  Distribution 
(See  figure  12  for  scaling) 


lgure  45a.  Flight  2-3-1  X  Axis  Probability  Density 
wee  Figure  12  for  scaling) 
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Figure  45b.  Flight  2-3-1  X  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 
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Figure  /*6b.  Flight  2-3-1  Y  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 


Figure  47a.  Flight  2-3-1  Z  4.>cia  Probability  Density 
(See  2'igure  12  for  scaling) 
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Figure  47b.  Flight  2-3-1  Z  Axis  Probability  Distribute 
(See  figure  i2  for  scaling) 
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Figure  A8b.  Flight  2-3-1  Roll  Acceleration  Probability  Distribution 
(See  figure  1 2  for  scaling) 
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Figure  49b.  Flight  2-3-.  X  /uti  ;  Probability  DjytribuMor 
(Jee  figure  12  for  scaling) 
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Figure  50b.  Fligh 
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Figure  52b.  Flight  2-3-2  Roll  Acceleration  Probability  Distribution 
(See  figure  12  for  scaling) 


Figure  S3a.  Flight  U 


re  •'  M  .  Flight  4-2  X  Axis  Probability  Distribution 
(Jee  figure  12  for  scaling) 


igure  54b.  Flight  4-2  Y  Axis  Probability  Distribution 
(See  figure  IP.  for  scaling) 


Figure  Flight  U-2  Z  Axis  Probability  Distribution 

(oee  figure  12  for  scaling) 


Figure  ;ba.  Flight  L-2  Roll  Acceleration  Probability  Density 
(£ee  figure  12  for  scaling) 


Figure  $6b.  Flight  4-2  Roll  Acceleration  Probability  Distribution 
(See  figure  12  for  scaling) 
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Figure  58b.  Flight  4-3  Y  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 


Figure  59a.  Flight  4-3  1  Axis  Probability  Density 
(See  figure  12  for  scaling) 


Figure  59b.  Flit  it  4-3  Z  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 
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Figure  60b.  Flight  4-3  Roll  Acceleration  Probability  Distribution 
(See  figure  12  for  scaling) 


Figure  6 If. .  Flight  U-U  X  Axis  Probability  Density 
(See  figure  12  for  scaling) 


Figure  61b.  Flight  4-4  X  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 


Figure  62a.  Flight  U-U  Y  Axis  Probability  Density 
(See  figure  12  for  scaling) 


Figure  63b.  Flight  4-4  Z  Axi3  Probability  Distribution 
(See  figure  12  for  scaling) 
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Figure  64b.  Flight  4-4  noli  Acceleration  Probability  Distribution 
(See  figure  12  for  scaling) 
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‘igure  65a.  Flight  6-2  X  4xis  Probability  Den3ity 
(See  figure  12  for  scaling) 
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Figure  65b.  Flight  6-2  X  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 


Figure  66a.  Flight  6 —2  Y  Axis  Probability  Density 
(See  figure  12  for  scaling) 
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Figure  66b.  Flight  6-2  Y  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 
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Figure  67a.  Flight  6-2  Z  Axis  Probability  Density 
(See  figure  12  for  scaling) 


Figure  6?b.  Flight  6-2  Z  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 
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Figure  68a.  Flight  6-2  Roll  Acceleration  hrcbability  Density 
(See  figure  12  for  scaling) 
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Figure  69a.  Flight  6-3  X  Axis  Probability  Density 
(See  figure  12  for  scaling) 


Figure  69b.  Flight  6-3  X  Axis  Probability  Distribntion 
(See  figure  12  for  scaling) 
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Figure  70b.  Flight  6-3  Y  Axis  Probability  Distribution 
(See  figure  12  for  scaling) 


Figure  ?lb.  Flight  6-3  Z  Axis  Probability  Distribute 
(See  figure  12  for  scaling) 
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Figure  72a.  Flight  e>-3  Roil  Acceleration  Probability  Density 
(See  figure  12  for  scaling) 
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Figure  74a.  Flight  1-1-2  X,  Y  and  Z  Axes  PSD 
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Figure  74b.  Flight  1-1-2  Roll  Axis  PSD 
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Figure  75R.  Flight  1-1-3  Roll  Axis  PSD 
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Figure  76b.  Flight  1-2-1  Ho 11  Axis  PSD 
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SWEPT  SPECTRUM  ANALYZED 
B=0 . 2  Hz,  T=3i*0  SEC. 


Figure  77b.  Flight  1-2-2  Roll  Axis  PSD 


Figure  78a.  Flight  1-2-3  X,  Y  and  Z  Aj^s  F3D 
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Figure  78b.  Flight  1-2-3  Roll  Axis  PSD 
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Figure  79b.  Flight  2-1  Roll  Axis  PSD 
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Figure  81b.  Flight  2-3-1  Roll  Axis  FSD 
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Figure  82a.  Flight  2-3-2  X,  Y  and  2  Axes 
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Figure  84b.  ’■'light  4-3  Roil  Axis  PSD 
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igure  &t> b.  Flight  4 -A  Holi  Axis  PSD 
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Figure  87a.  Flight  6-3  X,  Y  and  Z  Axes  PSD 


Figure  87b,  Flight  6-3  Roll  Axis  PSD 
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Figure  88.  Composite  X  Axis  PSD's 
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SPECTKUM  ANALYZED 


Figure  98.  Comparison  of  Experimental  -  Hypothetical  Y  Axis  PSD's 


Figure  99.  Comparison  of  Experimental  -  Hypothetic, 
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